Magnetic properties of a-iron(II) phthalocyanine 
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We report on the magnetic properties of the supra-molecular compound iron(II) phthalocyanine 
in its a-form. dc- and ac-susceptometry measurements and Mossbauer experiments show that the 
iron atoms are strongly magnetically coupled into ferromagnetic Ising chains with very weak antifer- 
romagnetic interchain coupling. The transition to 3D magnetic ordering below 10 K is hindered by 
the presence of impurities or other defects, by which the domain-wall arrangements along individual 
chains become gradually blocked/frozen, leading to a disordered 3D distribution of ferromagnetic 
chain segments. Below 5 K, field-cooled and zero-field-cooled magnetization measurements show 
strong irreversible behavior, attributed to pinning of the domain-walls by the randomly distributed 
defects in combination with the interchain coupling. High-field magnetization experiments reveal a 
canted arrangement of the moments in adjacent ferromagnetic chains. 

PACS numbers: 75.50.Xx; 76.80.+y 



I. INTRODUCTION 

Metal phthalocyanines (hereafter referred to as MPc) 
have received extensive attention in the last three decades 
because of their commercial applications as dyes, cat- 
alysts, and because of their interesting electro-optical 
properties [Q 0|. Their magnetic properties have at- 
tracted attention as well since manganese phthalocyanine 
(MnPc) was one of the first molecular magnets || . In this 
work we focus on iron(II) phthalocyanine (FePc). 

In the MPc molecule, the M atom has square planar 
coordination with four pyrrolic N atoms. In the bulk and 
in thin films, these planar molecules form linear chains 
by columnar stacking. The molecules are tilted by an 
angle with respect to the chain axis j|, ||. The angle at 
the parallel adjacent chain is opposite, thus constructing 
a so-called "herringbone" structure. The different angles 
and relative positions of the chains give rise to different 
polymorphs, for example in the a and /3 forms of copper 
phthalocyanine (CuPc) ||. This is of importance since 
the magnetic properties of metal complexes are strongly 
dependent on structural differences. For example, the 
particular structure in the /3-MnPc is thought to underly 
the Mn-Mn ferromagnetic interactions that make this 
compound the only ferromagnetic MPc compound so far 
reported ||. 

A very similar molecule with square planar M coordi- 
nation is M-porphyrine. These molecules may also stack 
in different polytypes, giving rise to different magnetic 
properties. Very recently it was reported that Fe(II) oc- 



taethyltetraazaporphyrin crystallizes in two structures: 
the a-form which is a canted, soft molecular ferromag- 
net with Tc — 5.6 K, and the /3-form which remains 
paramagnetic down to 1 K || Apparently, though 
the Fe(II) electronic states are identical, the different 
exchange paths and distances between the neighboring 
molecules and between chains modulate the magnetic be- 
havior. Thus, also Fe(II) square coordinated complexes 
seem prone to show various types of magnetic behavior. 

FePc could fall into this category of magnetic com- 
plexes. There are two forms of FePc described in the lit- 
erature. The a-form is metastable and obtained either as 
a polycrystalline powder or as a thin film formed by vac- 
uum deposition onto a cold substrate || B[. The /3-form 
is the most stable one and is obtained not only from the 
sublimation technique in the form of a single crystal, but 
also from heating the a-form sample to 350° C for a few 
hours Since 1935, the /3-form of this compound has 
been repeatedly reported to behave as a singlet-ground- 
state paramagnet down to the lowest measured tempera- 
tures [h], |ll], [l^] . On the other hand, to the best of our 
knowledge, the magnetic properties of a-FePc have so far 
never been studied. We show below that in the a-form 
of FePc the Fe(II) moments are strongly coupled into 
ferromagnetic chains, with the weak interchain coupling 
leading to a a canted, soft molecular ferromagnet below 
about 10 K. Our findings demonstrate that FePc has a 
behavior very similar to MnPc, so the magnetic behavior 
of the latter is not so exotic as thought earlier J3], but 
is a member of a more general category of square-planar 



Typeset by REVTeX 



2 



metal complexes in which ferromagnetic correlations are 
established. They represent, therefore, quite interesting 
examples of molecular magnets, a subject of intensive 
on-going research [fl3|| . 



II. EXPERIMENTAL DETAILS 

All the experiments were performed on powdered sam- 
ples of the a-FePc compound as obtained from Aldrich 
company (catalogue number 37,954.9), in two different 
batches. We verified the stoichiometry by elemental anal- 
ysis. The morphology was studied by scanning electron 
microscopy with a resolution of 3.5 nm. A step scanned 
powder X-ray diffraction pattern was collected at room 
temperature using a rotating anode X-ray diffractometer. 
The diffractometer was operated at 80 mA and 40 kV 
with Cu radiation and a graphite monochromator was 
used to select the Cu Kai,2 radiation. Data were col- 
lected from 4° up to 90° with a step size of 0.02° and a 
counting rate of 4 s/step. 

Magnetic data were obtained with the use of a commer- 
cial SQUID based ac-magnetometer instrument down to 
1.7 K with an external magnetic field in the range —50 
to 50 kOe. The amplitude of the ac-field was 4.5 Oe 
in all the experiments and the frequency / was varied 
between 0.1 and 1380 Hz. Magnetic data with larger 
applied field (up to 200 kOe) were obtained in the Ni- 
jmegen High Field Magnet Laboratory. The magneti- 
zation was measured using the ballistic extraction tech- 
nique and the sensitivity was about 10 -3 emu. The heat 
capacity was measured from 1.8 K up to 20 K with a 
commercial calorimeter. 



III. MORPHOLOGY AND STRUCTURE 

By scanning electron microscopy, we have observed 
that a-FePc develops in crystallites whose shape is sim- 
ilar to a 'bunch of rods' (Fig. 1). The ratio of length 
to thickness is approximately 10. Without any preferred 
length, the crystallites can be as long as up to hundreds 
of micrometers. By means of powder X-ray diffraction we 
have found that the present form of FePc is isomorphic to 
the al-form of 4-monochloro-copper phthalocyanine 
that is triclinic with space group PI or PI. The fit of the 
X-ray diffraction spectrum provides the following cell di- 
mensions: a = 12.2 A, b = 3.78 A, c = 13.0 A, a = 90.0°, 
(3 = 90.4° and 7 = 95.6°. The fit to the peak intensities 
is not good due to the needle shape of the powder grains 
that produces a texture effect. Thus, in the present sam- 
ple of a-FePc the molecules are stacked parallel to each 
other along the short &-axis at intervals of 3.78 A, and the 
molecular planes arc inclined to the ac plane at an angle 
of 26.5°, in a herringbone structure (Fig. 2). The verti- 
cal distance between the neighboring molecules is about 
3.4 A. We have verified further that the present form of 
FePc is indeed the (metastable) a-form because it shows 



the characteristic IR-absorption peak at 250 cm' 1 M. 
By warming our sample at 350° C during 24 hours in 
vacuum, the compound transforms into the /3-form, as 
expected. 

In the a-form, the Fe atom is coordinated to the four 
N atoms in the diagonal positions of the square pla- 
nar molecule. The nearest neighbors on the adjacent 
molecule at one side are the two isoindole and one aza- 
methine nitrogens located at similar distances, while at 
the other side the coordination is identical but rotated by 
180°. The Fe(II) ion is in a symmetry lower than D^, 
but in a first approximation we shall assume this local 
symmetry. The value of 3.78 A for the Fe-Fe separation 
along the 6-axis in the present form of FePc is very rea- 
sonable, lying within the accepted range of 3.5-4.8 A for 
metalloporphyrin dimers fl4|| . It is generally accepted 
that the Fe-Fe separation along the 6-axis and the an- 
gle of inclination of the molecular planes with respect to 
the crystallographic ac plane determine the extent of the 
exchange Fe-Fe interaction. 

We note that the main structural difference with the 
/3-form consists in the inclination angle (Fig. 2); in the 
/3-form it is 45.8° and the azemitine nitrogen of each ad- 
jacent molecule lays just on the apical position of an oc- 
tahedron of N atoms ||, in the a-form it is 26.5°, and 
the N atoms from the nearest molecules are not in axial 
positions and therefore do not form such a N octahedron 
about the Fe atom. 



IV. EXPERIMENTAL RESULTS 

The inverse of the in-phase component of the zero-field 
ac-susceptibility of a-FePc, as measured at / = 90 Hz, is 
plotted as a function of temperature in Fig. 3. As shown 
by the solid line, the data can be fitted between 50 K and 
300 K by a Curie- Weiss law X = C/(T - 6) for 5 = 1, 
g = 2.54 and = (40 ± 2) K. This result implies that in 
this temperature range the iron(II) ion is in an effective 
spin 5 = 1 state, and the large positive value of 9 in- 
dicates that there are strong ferromagnetic correlations 
between the Fe moments in the paramagnetic phase. On 
basis of the powder X-ray structural analysis, we expect 
these correlations to extend within the columnar chains 
along the &-axis. Neglecting interchain magnetic coupling 
and using the mean-field expression for the Curie- Weiss 
temperature 6 = 2zJS(S+l)/3ks, we find for 5 = 1 and 
z = 2 the estimate J/k-Q m 15 K for the ferromagnetic 
interaction strength along the chains. We note that also 
for the /3-form of FePc, an effective spin 5=1 was found 
below room temperature. Indeed, at high temperatures, 
the susceptibilities of the a and j3 ||| [l5| forms are coinci- 
dent (inset of Fig. 3). The fact that the susceptibility of 
the /3-form goes to zero for T — » has been explained in 
terms of zero-field splitting of the triplet state, leaving a 
non-magnetic singlet as the lowest energy level occupied 
for T< 20 Kg. 

Turning now to the data taken below 50 K, as exempli- 
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fled in Fig. 4, we note that the behavior of the susceptibil- 
ity of a-FePc points to strong low-dimensional magnetic 
characteristics |L6|. Indeed, as seen in Fig. 4, the suscep- 
tibility continues to increase down to temperatures much 
lower than 8. Only at about 10 K the ac-susceptibility 
curve shows a change of slope, and a rounded maximum 
is found at still lower temperatures. The paramagnetic 
behavior is thus sustained down to approximately 10 K, 
that is much less than 9, as expected for a magnetic chain 
system Jl6| . 

From the shape of the X'(^) curve in Fig. 4, one would 
be tempted to associate the change of slope to the onset 
of 3D correlations between the chains due to antiferro- 
magnetic interchain coupling. However, here we should 
point to the fact that \' at the maximum is not very far 
below the calculated limit for a completely ferromagnet- 
ically ordered isotropic material below Tc, namely 1/N 
with N the demagnetizing factor (a rough estimate gives 
the value of w 33 emu/mol). 

Since the presence of anisotropy can lower the value 
of the susceptibility of a ferromagnetic powder consider- 
ably Jl6| , |l7j ], the occurrence of ferromagnetic interchain 
ordering would also be considered. Moreover, in case the 
moments of the ferromagnetic chains would be subject to 
strong crystal field anisotropy, the latter could compete 
with the interchain interactions, and, when sufficiently 
strong, lead to a freezing at some blocking temperature 
before the interchain coupling becomes effective (analo- 
gous to the phenomena observed in molecular magnetic 
clusters and superparamagnetic particles). 

In order to further investigate the nature of the "order- 
ing" processes, we have therefore measured the zero-field 
differential susceptibility in the region below 15 K for 
different frequencies. The dependence on frequency of 
the real and imaginary parts of the ac-susceptibility are 
shown in Fig. 5. As was already seen for / = 90 Hz in 
Fig. 4, for all frequencies used the in-phase component 
of the ac-susceptibility, x', shows a rounded maximum 
in the range 5 K-10 K. This is always accompanied by a 
maximum in the out-of-phase component, x", at some- 
what lower temperature. The maximum in x'(T) shifts to 
lower temperature with decreasing frequency, and x'(T) 
decreases abruptly to zero below the maximum. The 
sharp maximum in x"{T) shifts likewise with frequency. 
It is not possible then to define a definite transition tem- 
perature to a 3D ordered state from these data, since the 
obtained Tc would be frequency dependent. The shift in 
temperature of the maxima in x'(T) and x"(T) with the 
frequency are reminiscent of spin-glass or superparam- 
agnetic behavior. Since the freezing temperature Tf(u>) 
can be defined from wT[Tf(w), H = 0] = 1, a quanti- 
tative measure of the frequency shift is obtained from 
(AT f (uj)/T f (uj))/A\n (u) = 0.034. It is six times larger 
than the rate for a metallic spin glass and an order of 
magnitude smaller than for a superparamagnet Jl8| . 

In view of the just-described frequency dependence of 
the ac-susceptibility, we decided to investigate the anal- 
ogy with spin glass systems or superparamagnets still 



further, by means of dc-magnetization studies. In Fig. 6 
we show dc-M(T) measurements taken with the SQUID 
magnetometer in the zero-field cooling (ZFC) and field- 
cooling (FC) modes at two fields, H = 50 Oe and 1 kOc. 
We note that both curves show a change in slope at about 
10 K, similar as for x' hi Figs. 4 and 5. Furthermore, a 
bifurcation point showing strong irreversibility is found 
below about 5 K. At first sight, the bifurcation tempera- 
ture is not much affected by the field although we point 
out that for the low-field curve (50 Oe) the ZFC and FC 
curves strictly speaking only coincide above ~ 10 K. 

In another experiment, after a FC process down to 
1.8 K with H = 40 Oe the field was switched off and 
the remnant magnetization was measured (Fig. 7). We 
found a large remnant magnetization at 1.8 K which very 
rapidly decays by several orders of magnitude when in- 
creasing the temperature towards 5 K. Above 5 K, the 
decay of the remanence proceeds much slower (inset of 
Fig. 7), and a small amount of remnant magnetization is 
seen to persist even up to 10 K-15 K. On basis of Fig. 7 
we may conclude that, although most of the field- induced 
magnetization is indeed lost at about 5 K, a small frac- 
tion of about 3% is left above this temperature, and is 
only removed by warming up to about 15 K. This behav- 
ior points to the presence of substantial ferromagnetic 
short-range order along the chains in this intermediate 
region of 5 K-15 K, in agreement with the behavior of 
the ac-susceptibility. 

The M{H) measurements of the initial (ZFC) mag- 
netization were therefore extended to higher fields and 
at temperatures above, at and below this intermediate 
region in order to get some more insight in these irre- 
versibility processes (Fig. 8). The data at T ~ 30 K 
show the expected paramagnetic behavior. In the in- 
termediate temperature region, for T = 8 K and 5 K, 
the initial increase in M{H) is very abrupt, reaching the 
value 1.4 \xb for H = 2.5 kOe. Above 10 kOe, the magne- 
tization curves measured with the SQUID up to 50 kOe 
(Fig. 8) or in the Nijmegen High Field Magnet Labo- 
ratory (Fig. 9), show an additional slow increase up to 
a = 1.7 [i B at H = 200 kOe. Note in Fig. 9 the per- 
fect coincidence of the SQUID equilibrium measurement 
with the branch corresponding to the high-field measure- 
ment at decreasing field, while the branch measured at 
increasing field is lower. The SQUID measurements have 
a characteristic point-to-point measuring time of ~ 150 s 
while the high-field measurements have a characteristic 
time of ~ 40 s. The system evidently needs a longer time 
than 40 s to reach equilibrium at this temperature when 
the field is increased from zero. 

The spontaneous magnetization evaluated as the ex- 
trapolation to H = of the high field demagnetizing 
curve yields a value of 1.5 \xb- We point out that, al- 
though these measurements demonstrate that in the re- 
gion below about 5 K the compound is in a ferromagnetic 
state, when the g factor deduced in the paramagnetic re- 
gion is still valid in this range, the saturation magneti- 
zation would be 2.5 /ij on basis of the (effective) spin 
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5 = 1 ground state found for the Fe(II) ions from the 
Curie- Weiss analysis of the high-temperature suscepti- 
bility (Fig. 3). We shall come back to this discrepancy 
below in the discussion section, where we shall argue that 
on basis of the competition between the exchange term 
and the crystal field term in the Hamiltonian for a-FcPc, 
a substantial canting of the magnetic moments with re- 
spect to the crystal field (anisotropy) axis is to be ex- 
pected. The increases seen in the M(H) curves above 
5 kOe may then be associated with the gradual rotation 
of the canted moments, produced by the competition be- 
tween the Zeeman energy and the crystal field energy. 

One may conclude that below 5 K, the material shows 
strong irreversibility, as already evidenced by the bifur- 
cation point in Fig. 6. Also the behavior of the M(H) 
curves in Fig. 8 points to the presence of long relaxation 
times. For example, the initial M(H) curve measured at 
1.8 K is actually even lower than the data at 15 K, since 
the time waited to measure each point is of the order 
of minutes and only a fraction of the sample is oriented. 
This is more clearly seen in the hysteresis loop measured 
at 1.8 K by sweeping the field between —50 and 50 kOe 
in a time of ~ 150 s (Fig. 10). The initial curve increases 
progressively with the field, tending to a value of 1.6 \xb 
at 20 kOe. When the field decreases, the magnetization 
remains nearly constant down to zero-field, yielding a 
remanence of 1.5 [Ib- However, the inversion of the rem- 
nant moment is very abrupt when a field in the opposite 
direction is applied, indicating that there is no coercivity 
in the usual sense associated with this remanence. 

Similarly strong relaxation effects in the low- 
temperature range have been seen in preliminary 57 Fe 
Mossbauer spectroscopic studies on a-FePc. In Fig. 11 
spectra at some representative temperatures are shown. 
Above 25 K, the spectrum is hardly temperature depen- 
dent and reveals the familiar paramagnetic, quadrupolar 
split doublet. Below that temperature the doublet starts 
gradually to broaden and below about 15 K a magnetic 
hyperfine splitting appears, at first extremely broadened, 
then becoming the sharper the lower the temperature. 
Below about 5 K, the broadening has almost disappeared 
and a fully hyperfine split spectrum is observed. A de- 
tailed presentation of the Mossbauer spectra will be given 
in another publication. 

Similar temperature dependent relaxation phenomena 
as seen in the Mossbauer spectra in Fig. 11, have typically 
been seen for blocking phenomena in superparamagnetic 
particles or for domain-wall excitations in anisotropic 
magnetic (Ising) chains. Briefly, the magnetic moment 
of the Mossbauer nucleus (here 57 Fe) feels the hyperfine 
field exerted on it by the electron spin on the same atom. 
Thus the fluctuations in the hyperfine field are directly 
related to the electron spin dynamics, i.e. the nuclear 
moment probes the autocorrelation function of the elec- 
tron spin. Accordingly, the Mossbauer spectrum will be 
broadened as long as the fluctuation rate of the electron 
spin lies in the frequency window (10 6 Hz-10 9 Hz) of 
the Mossbauer experiment. The temperature dependence 



observed indicates that the electron spin dynamics is de- 
scribed by thermally activated processes. At high tem- 
peratures the fluctuation rate of the electron spins is too 
fast to cause broadening. With lowering temperature the 
rate slows down, enters the Mossbauer frequency window 
leading to the strong relaxation. At the lower tempera- 
tures the rate becomes just too slow to be felt. In super- 
paramagnetic particles, the activated process is related 
to the thermally excited switching of the particle's mag- 
netic moment between the orientations dictated by the 
anisotropy, which becomes blocked at low temperatures. 
In the magnetic Ising chain systems, the fluctuation rate 
is dictated by the density of thermally excited domain- 
walls (solitons) along the chains. At high temperatures, 
the density is very high and thus the fluctuations too 
fast. At intermediate temperatures, the fluctuation rate 
passes the Mossbauer window, whereas, at low temper- 
atures, the soliton density becomes too small, moreover 
the walls may become blocked/pinned by interchain in- 
teractions or defects in the chains. 



V. DISCUSSION 

Summing up the evidence obtained so far, we may state 
that the magnetic structure of a-FePc is most probably 
composed of very weakly coupled ferromagnetic chains, 
with an intrachain exchange J/fce ~ 15 K, and an ef- 
fective spin 5 = 1 for the Fe(II) ion in the temperature 
range 50 K-300 K. The crystal field splitting, responsi- 
ble for the effective spin 5 = 1, likely result in addition 
into a strong anisotropy of the ferromagnetic interaction. 
Below, we shall thus start the discussion with a full anal- 
ysis of the crystal field effects in this material, in partic- 
ular also in relation to the other members of the metal 
phthalocyaninc and metal porphyrine family. Next we 
shall present the ensuing analysis of the magnetic suscep- 
tibility in the range 10 K-50 K in terms of ID magnetic 
models systems. Finally, the nature of the relaxation and 
ordering phenomena observed below 15 K will be further 
analyzed. 

As is well known, for Fe(II) the combination of or- 
thorhombic crystal field splitting combined with spin- 
orbit coupling leads to an effective spin description for the 
five lowest occupied levels (d xy , d yz , d xz , d z 2, d x 2_ y 2) |ljj, 
pp[ . Depending on the ratio of the parameters, and the 
relevant temperatures of the experiment, one can have a 
situation where either two, three or all five levels will par- 
ticipate in the magnetic ordering processes, leading to an 
effective spin 5 = 1/2, 1 or 2, respectively, for the Fe(II) 
ion. For thermal energies in the range 100 K-300 K, the 
a-FePc may be described as an effective spin 5 = 1 with 
g = 2.54. However, the thermally occupied levels, which 
lead to the 5=1, are split into a lower- lying singlet 
and an excited doublet, with the splitting being probably 
of the order of 50 K. Accordingly, the magnetic behav- 
ior will be determined by the competition between the 
magnetic exchange energy and the crystal field term rep- 
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resenting this singlet-doublet splitting. From the Curie- 
Weiss fit to the susceptibility, we have seen that 9w40K, 
so both these terms are of the same order of magnitude. 
Below we will show that the competition may actually 
produce a canted arrangement of the ferromagnetically 
aligned spins, with an effective spin S — 1/2 at low tem- 
peratures and considerable anisotropy. 

For the temperature range 100 K-300 K, we con- 
clude that the Fe(II) ion in a-FePc may be described 
as an effective spin S — 1 compound, in which single-ion 
anisotropy and second-order spin-orbit coupling leaves 
the triplet as the lowest energy multiplet. This is the 
same spin state as found in /3-FePc j|. Indeed, both 
compounds have identical magnetic properties at room 
temperature, which arise from a similarly split triplet 
[-D/&B — 53 K and 98 K, for a and /3-form, respectively 
(see below)] with the non-magnetic singlet as the ground 
state |pf. 

In the a-form the S" = 1 spins remain ferromagnetically 
coupled in magnetic chains down to the lowest temper- 
atures in spite of the non-magnetic ground singlet since 
the ferromagnetic interaction 2zJ (with J//cb ~ 15 K 
as estimated from the Curie- Weiss fit) is larger than D, 
in agreement with Moriya's criterion for magnetic order- 
ing in a S = 1 triplet pifl . Apparently, this condition is 
not satisfied in /3-FePc because the ration of D over J is 
larger, and as a consequence it remains paramagnetic at 
all temperatures. 

Curiously, the magnetic phenomenology of the MnPc 
forms is just opposite: i.e. the /3-form is ferromagnetic 
while the a-form is antiferromagnetic p2| . In the Mn 
compounds the different behavior of the two forms was 
attributed to competing ferro- and antiferromagnetic ex- 
change interactions, neglecting direct Mn-Mn coupling. 
Under the assumption that the electronic configuration 
of the 3d 5 electrons is (d xy ) 2 (dxz) 1 (dyz) 1 (d^) 1 , the ferro- 
magnetic interaction is caused by the overlap between the 
d z 2 unpaired electron (a\ g symmetry) with the 7r orbital 
of Eg symmetry of the adjacent aromatic ring, which is 
orthogonal to the d z 2 orbital of its metal ion (Fig. 12). 
The antiferromagnetic coupling is due to the overlap of 
the d xz d yz orbitals (e g ) with the same 7r orbital of E g 
symmetry (Fig. 12). In the MnPc /3-form the ferromag- 
netic mechanism dominates while in the MnPc a-form 
the ferromagnetic path weakens and antiferromagnetism 
prevails ||^]. 

If we assume that the electronic configuration of Fe is 
related to that of the Mn by just adding one electron 
to the second level as proposed for the /3-form of FePc 
by Dale et al. |l2j and Coopens et al. |^3j, that is E g A 
{d X y) 2 (d X z) 2 (dyz) (d z 2) 1 , then the same arguments as for 
the Mn compounds would be applicable. The only dif- 
ference would be that the antiferromagnetic interaction 
should be weakened, since the (d xz ,d yz ) doublet should 
have just one unpaired electron. Consequently, the /3- 
form should be ferromagnetic, contrary to observation. 
One arrives to the same conclusion if one assumes the 
electronic configuration £?2 S (dxy) 1 (d X z) 2 (dyz) 2 (dz 2 ) 1 , as 



proposed by Barraclough et al. Q. 

The conclusion is that the electronic configuration 
needs to be different, or direct exchange is more intense 
in the a-form than in the /3-form. It can be argued that 
the orbital ground state A 2g A (dxy) 2 (dxz) 1 (d yz ) 1 (d z 2) 1 , 
proposed by Stillman et al. (24| leads to just antiferro- 
magnetic coupling. The remaining configuration E g B 
(d z 2) 2 (d xz ) 2 (dy Z ) 1 (d X y) 1 , proposed by Labarta 0, gen- 
erates a path for ferromagnetism in the a-form: the un- 
paired electron in the d xy orbital singlet may interact 
with its counterpart in the adjacent molecule via the E g 
it electron of the aromatic ring, which is orthogonal, giv- 
ing rise to its ferromagnetic character (Fig. 12). This is 
the only ferromagnetic exchange path possible, however 
its intensity would be probably very weak since it is in 
the xy plane and the overlap will be small. At any rate, 
in the /3-form the d xy orbital has negligible overlap with 
the E g electron of the apical N since it is located in the 
z axis, thus the ferromagnetic interaction would be neg- 
ligible. The unpaired electron of the (d xz ,d yz ) doublet 
would lead in both forms to antiferromagnetic interac- 
tion. So, exchange interaction is not likely to cause the 
ferromagnetic coupling present in the a-form. 

On the other hand the condition for an enhancement 
with respect to the /3-form of the direct Fe-Fe inter- 
action are fulfilled since the Fe-Fe distance is reduced 
from 4.8 A(/3) to 3.8 A(a). Following Goodenough ||, 
the direct ferromagnetic interaction occurs when there 
is overlap between a half filled and an empty orbital, 
or a half-filled and a full orbital. Since the (d x 2_ y 2) 
empty orbital is very high in energy |l5| we think that 
only the case of half-filled full-orbital overlap may be rel- 
evant here. Then, if the orbital configuration is E g A 
(d xy ) 2 (d xz ) 2 (dyz) 1 (d^) 1 , as proposed by Coppens et al. 
for the /3-form from accurate X-ray diffraction |^3), the 
overlap may take place between two (d xz ) 2 (d yz ) 1 orbitals. 
The proposition of the E g A electronic configuration by 
these authors is based on very direct observation of elec- 
tronic deformation densities and in turn supports the 
proposed mechanism for the ferromagnetism in the a- 
form. Of course, this overlap will decay strongly when 
the atoms are separated further by 1 A, and so the fer- 
romagnetic interaction weakens to the extent of being 
negligible in the /3-form. 

Thus, we conclude that the difference in the phe- 
nomenology of the Mn and Fe forms can be related to 
the different ground state spin configuration (two mag- 
netic doublets for Mn, and one magnetic doublet and a 
non-magnetic singlet for Fe) on the one hand, and to the 
different electron filling scheme and overlap probabilities 
of the resulting orbitals, on the other. 

As said in the introduction, such a different behavior, 
ferromagnetism or paramagnetism, between two herring- 
bone structured compounds formed by stacked square 
planar molecules containing Fe, has also been detected 
in Fe(II) octaethyltetraazaporphyrin || Similarwise, 
the a-form is ferromagnetic while the /3-form is paramag- 
netic. However, in a recent paper |2q] a magnetic study 
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of Mn(II) octaethyltetraazaporphyrin shows that the a- 
form is also ferromagnetic while the /3-form is paramag- 
netic, in contrast to the behavior of the MnPc forms. 
So, the parallelism between the phthalocyanines and oc- 
taethyltetraazaporphyrin is not complete. 

On the basis of the above, it appears to be beyond 
doubt that the magnetic moment of the Fe(II) in the 
a-FePc is to be associated with a spin triplet probably 
split by crystal field into a doublet and a lower-lying sin- 
glet. We then have to reconcile this conclusion with the 
rather low values of the magnetization reached even in 
high fields as exemplified in Figs. 8-10. Based on an ef- 
fective spin S = 1 and the powder g- value of 2.54 deduced 
from the Curie- Weiss fit to the susceptibility above 50 K, 
the saturation moment should amount to 2.5 /xb per ion, 
i.e. far above the measured value even for an applied 
field of 200 kOe. Looking more closely at the measured 
M(H) curves and neglecting for the moment the relax- 
ation effects, one may observe that below about 10 K the 
magnetization process can be divided in two steps. From 
the M(H) curve measured at 5 K (Fig. 8), a very rapid 
increase of the magnetization up to 1.0 [Ab (per atom) 
is observed for fields up to 400 Oe, followed by a much 
slower gradual increase up to the highest field applied. 
This two-step behavior can be explained in terms of a 
division of the magnetic chains into two different mag- 
netic sublattices, with a mutually canted arrangement. 
In that case the fast rapid increase would correspond to 
saturation of the moments within each of the magnetic 
sublattices, followed by the gradual rotation of the sub- 
lattices to a common alignment parallel to the field. In 
the latter process, the Zeeman energy has to compete 
with the (large) crystal field anisotropy energy, so that 
the high-field susceptibility (slope of the M(H)) curve 
can be very small. It is of interest to point out in this 
regard that for /3-MnPc a canted magnetic structure in 
the above sense has indeed been reported [^7| . 

At first sight the preposition that the ordered magnetic 
structure is canted seems to contradict the positive value 
of D, which implies a single crystal anisotropy of planar 
character. Indeed, Miyoshi et al. |^7| concluded that for 
chains consisting of ferromagnetic coupled Fe moments 
with isotropic exchange, and with single-ion anisotropy, 
such that the quantization axis of the Fe sites in one chain 
is tilted by an angle which is opposite respect to the tilt 
angle in the adjacent chain, a negative D is needed to 
stabilize a canted structure, irrespectively if the inter- 
chain interaction is ferro or antiferromagnetic. This pre- 
diction, when applied to the MnPc compound, revealed 
a contradiction between the canted structure found in 
that compound by magnetization measurements and po- 
larized neutron diffraction H , and the positive D deter- 
mined earlier from the analysis of the high-temperature 
susceptibility a discrepancy that has gone unex- 

plained so far. This apparently identical contradiction, 
now found for a-FePc and /3-MnPc, can be explained 
by noting that, contrary to the assumption of Miyoshi 
et al. p7[ that the ferromagnetic intrachain interaction 



is isotropic, we find from our experiments that it is in- 
stead of the anisotropic Ising-type, with the easy axis 
perpendicular to the molecule plane. Interestingly, the 
same result; i.e. Ising intrachain interaction, was found 
to fit best the MnPc high temperature susceptibility by 
Barraclough et al. [^8|, although it was not exploited. 
Thus, the Ising character of the interaction and its larger 
strength, comparable to the crystal field term, creates an 
uniaxial exchange anisotropy that opposes the crystal- 
field anisotropy. A rough mean field calculation for FePc, 
using the value J/fce ~ 15 K found from the Curie- Weiss 
9, tells us that H ex = (2z,J/k B ){S z )/gfi B = 58 T (at 
T = K). Such an intense field splits the upper S z = ±1 
doublet and brings the S z = — 1 level 30 K down so that 
it nearly coincides with the S z — level. The net effect 
is that, in terms of the 5* = 1 formalism, the magnetic 
moment is polarized towards the quantization direction 
(local z axis) at low temperatures, which is the neces- 
sary condition to have the canted structure according 
to Miyoshi et al. |27j. In MnPc the exchange field is 
H ex = 48 T (D/k B = 29 K), and as a consequence the 
S z = —3/2 state is pulled below the S z = —1/2 state by 
68 K (at T = K) and the ground state gives rise to an 
effective uniaxial anisotropy. Thus the annoying contra- 
diction of an apparent inversion of the S z = —3/2 and 
S z = —1/2 states above and below the transition, men- 
tioned at the end of Ref. ||, is resolved. In the paper 
by Reiff et al. [^6| on Fe-porphyrine, a similar qualitative 
argument of opposing exchange and crystal field effects 
is mentioned to explain the magnetic properties of those 
compounds. 

The logical conclusion from the model we propose is 
that each chain has its own ferromagnetic axis tilted from 
the fo-axis by a given angle, with the easy axis in the 
adjacent chain tilted by the opposite angle, for reasons 
of symmetry. Thus, in a domain formed by this ordered 
magnetic structure there is a direction of uncompensated 
moments, and a perpendicular direction of compensated 
moments. The uncompensated moments should have the 
direction of the 6-axis in case the interchain interaction 
is ferromagnetic, and perpendicular to b if it is antiferro- 
magnetic. In either cases, there are two degenerate types 
of domains, in contrast to the four of MnPc, due to the 
different symmetry. We refrain from pushing any fur- 
ther the analysis since we have insufficient experimental 
information on the magnitude and sign of the interchain 
interaction nor on the precise tilt angle, because our mea- 
surements are performed on powder samples. 

We remark that the above analysis points to two dif- 
ferent scenarios for analyzing the magnetic susceptibility 
data below 50 K, namely in terms of ferromagnetic Ising 
chains with either effective spin S — 1, or effective spin 
S = 1/2, when considering only the lowest two levels 
(S z = and S z = —1). Below we shall investigate the 
applicability of both these possibilities. 

Starting with S = 1, we have calculated the parallel 
susceptibility of an S = 1 linear chain under the assump- 
tion of an Ising-like exchange interaction and of an uni- 
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axial single-ion anisotropy D as described by the Hamil- 
tonian: 

N—l N N 

Ha = -2 Ji ^ Sz,iS z ,i+i + D 2^ 5 2 ^ - g^ B H z 2J 

i=l i=l i=l 

The eigenvalues A±, Ao of this Hamiltonian were ob- 
tained by the conventional transfer matrix technique, 
similarly as described in Ref. |29|j , yielding: 

A± = \{l + a + 1) ± i[(a + 7 - I) 2 + 8/3 2 ] 1 / 2 (5.1) 

A = a - 7, 

where a = exp[(2J x - D)/k B T], 7 = exp[(-2J x + 
D)/k B T] and /3 = exp(D/2fc B T), with A± being the 
largest eigenvalue. Utilizing a Taylor series expansion 
to obtain the eigenvalues in the non-zero-field case, the 
susceptibility can be expressed as 

jVg 2 /i|^+l)e 2 VA + -l + e^ 

x = p 5 - 2 ) 

where 

F A = 3A 2 h -2A + (l + 2e-''+ 2d )-[2e 2d (l-e''V2e 4 ' i sinh(2j)], 

(5.3) 

with j = 2J l /k B T and d = -D/2k B T. 

Least squares fits of Eq. (^^) to the experimental sus- 
ceptibility data were carried out. Since the contribution 
of the perpendicular susceptibility has been neglected, 
only the data for relatively high temperatures (T > 10 K) 
were used. The dotted line in Fig. 4 shows the best fit 
with Ji/fc B = 25.7 K, g = 2.54 and D/k B = 53.2 K, 
where the positive J\ indicates ferromagnetic interaction 
and the positive D implies that the S z = ±1 doublet lies 
above the S z = state. The results are sensitive to both 
the temperature range utilized and the g value assumed. 
We have tried to fit the same data with the isotropic 
5 = 1 Heisenberg model and uniaxial anisotropy p0[ , 
but the fit is much worse (Fig. 4). For the temperature 
range from 20 to 300 K, the following result was obtained: 
J l /k B = 16.2 K and D/k B = 76.7 K. Since we are fit- 
ting powder data while the magnetic properties are very 
probably highly anisotropic, the values obtained from the 
fit have to be considered as approximate. However, we 
note that the anisotropy constant is about half the value 
of D/k B = 98 K found in the /3-form fl§]. Obviously, the 
different anisotropy originates from the absence of the 
apical N atoms next nearest to the Fe in the a-form re- 
spect to the /3-form. The most important result, though, 
is that the Fe-Fe intrachain interaction is unambiguously 
ferromagnetic. The value obtained for J\ is in fact similar 
to that following from the Curie- Weiss fit. 

We next turn to an analysis of the paramagnetic sus- 
ceptibility below 50 K in terms of the effective spin 
5 = 1/2 formalism (Fig. 4), which would apply when 



either the action of a crystal field of orthorhombic or 
lower symmetry, or the combined action of crystal field 
and magnetic exchange, would produce an effective dou- 
blet as the only occupied energy levels at low tempera- 
tures. In that case, since the effective anisotropy for the 
doublet will be of the Ising-type, we may compare the 
experimental powder susceptibility data with the theo- 
retical predictions for the parallel (xii ) and perpendicular 
(x_l) susceptibilities for the ferromagnetic 5 = 1/2 Ising 
chain @: 

X|| = (N g y B /2J 1/2 )Kexp(2K) (5.4) 

X± = (A^iMlM/2)[tanh(X) + Kcosh- 2 (K)} (5.5) 

where K = Ji/ 2 /2k B T, from which we have calculated 
the powder susceptibility \p — (X|| + 2%_l)/3. We can 
consider also that due to the large expected anisotropy 
one has g\\ 3> g±, and that the saturation magnetization 
at low temperature is M p = ^(g\\ + 2g±)S/i B , for which 
the experiment yielded l.5/i B . Fitting both quantities 
(for T > 10 K) leads to the estimates Ji/ 2 /k B = 32.2 K, 
gii = 8.6 and g± — 0.2. The fit of the magnetic suscepti- 
bility data with the prediction from the above equations 
is shown in Fig. 4. It is difficult to judge whether the fit 
on the basis of the 5=1 model is better than that with 
the 5=1/2 prediction since: (a) the deviations seen for 
T < 10 K from the 5=1/2 curve may well be due to 
the interchain interactions, which are not taken into ac- 
count; (b) for the same reason the seemingly better fit 
with the 5=1 model (Fig. 4) may be fortuitous since 
it also neglects interchain interactions; (c) as a further 
complication the x± contribution had to be neglected in 
the 5 = 1 calculation since no theoretical prediction is 
available. 

Proceeding next to a discussion of the order- 
ing/blocking phenomena observed below 10 K, we can 
start from the assertion that we may describe a-FePc as a 
system of very weakly coupled ferromagnetic Ising chains 
with either effective spin 1/2 or 1. For such systems, 
it is well-known that the ID magnetic behavior is dic- 
tated by the excitation of domain-walls (solitons) along 
the chains. For strongly anisotropic chains, the width of 
the domain- wall becomes smaller than a lattice constant 
and we speak of "kinks" separating spin-up and spin- 
down regions. At least two types of excitations in such 
Ising chains are possible. Namely, the single-kink exci- 
tations at the ends of the chains (with excitation energy 
J) and the kink-pair excitations (with energy 2 J). The 
correlation length along the chain (and thus the magni- 
tude of the magnetic susceptibility) is determined by the 
number of kinks excited (£ oc oc exp (—J/k B T)). 

Thus we think the present system falls in the same class 
of Ising-type chains studied previously, a.o. by magnetic 
susceptibility and Mossbauer effect studies, in connection 
with the soliton problem (see e.g. Refs. (||, [jj, |35|, 
|37| , |38| , |39||). In fact, the behavior observed here appears 
to be analogous to that of FeCl2Py2, one of the ferro- 
magnetic Ising chains studied in these previous works. 
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However, for a-FePc, the interchain coupling turns out 
to be much weaker than in FeCl2Py2 and related com- 
pounds. 

As already mentioned when presenting the Mossbauer 
data in Fig. 11, the extremely strong broadening of the 
spectra observed in the range 5 K-20 K may indeed be 
explained in terms of a spin dynamics based upon the ex- 
citation of domain- walls along the ferromagnetic chains. 
Since the Mossbauer effect is sensitive to high-frequency 
fluctuations, the relaxation effects already start to have 
an influence at about 20 K, whereas relaxation effects in 
the low-frequency susceptibility appear only below 10 K. 
In principle, both the high- and low-frequency spin dy- 
namics observed can be related to a gradual slowing down 
of the spin fluctuations with temperature due to the de- 
crease of the density and the mobility of the walls, fol- 
lowed by (almost) complete pinning below 5 K. It should 
be realized that, even when pinned, small excursions 
around equilibrium positions of the walls may still re- 
main possible, which explains the possibility of simulta- 
neous observation of both the high- and low-frequency 
spin dynamics in the range between 10 K and 5 K. 

At this point it is important to recall the similarity 
of the here-observed relaxation phenomena with those 
known for superparamagnetic particles pointed out sev- 
eral times in the above. Indeed, one might envisage 
the ferromagnetic chains in a-FePc as superparamag- 
netic particles, switching their magnetization among the 
up and down directions along the local anisotropy axis. 
However, it should then be immediately realized that 
the rotation of the spins along a given chain cannot 
occur simultaneously, in view of the extremely large 
anisotropy barrier involved in such a process (i.e. the 
anisotropy per spin times the number of spins in a cor- 
related chain). In full analogy with the situation in fer- 
romagnetic nanowires p0[ , the rotation has to occur via 
the excitation of a small reversed domain, either at the 
chain end by a 7r-wall, or along the chain by a kink- 
antikink pair, in both cases followed by the propagation 
along the chain of the wall(s). It is clear then, that the 
mobility of the domain-wall(s) along the chains may be- 
come an important factor in determining the speed of 
the switching process. In this respect the magnetic inter- 
actions between the chains, which tend to establish full 
3D order between correlated chain segments as the tem- 
perature is lowered, can be expected to slow down the 
wall motions along the individual chains. It is tempting, 
therefore, to associate the slowing down of the fluctua- 
tions in the range between 10 K and 5 K, seen in both 
the Mossbauer spectra and the ac-susceptibility, as a re- 
flection of the gradual establishment of 3D magnetic or- 
der between the chains. Since no clear transition at a 
definite temperature Tq is indicated, the resulting "or- 
der" is apparently not long-range but with a substantial 
amount of randomness, as in a spin glass. Judging from 
the behavior seen in e.g. Figs. 6 and 7 of the low-field 
magnetization, this ordering process is about completed 
below 5 K. The source of randomness is probably a ran- 



dom pinning of the domain- walls on impurities along the 
chain, whereby the 3D magnetic correlations between the 
chains are easily broken up. In view of the instability of 
a-FePc in favor of the /3-FePc form, the occurrence of 
such impurities is not unlikely. 

The above scenario can also explain the fact that in our 
specific heat data measured in the range 1.8 K to 20 K, 
no sign of a singularity (as associated with 3D long-range 
magnetic ordering) was detected. The measured specific 
heat curve just shows a smooth continuous behavior as 
expected for a combination of phonon contributions and 
ID magnetic fluctuations within the chains. 

Adopting therefore tentatively the value of 5 K as rep- 
resentative temperature for the ordering process between 
the chains, we may estimate the interchain coupling J' 
either from the Oguchi's formula |j4l| , or from the mean- 
field argument given by Villain p2j| , in both cases ob- 
taining J'/ks « (10~ 3 - 10~ 2 ) K. This estimate may be 
confronted with the behavior of the magnetization curves 
taken at 5 K and 8 K (Fig. 8), in which as discussed be- 
fore a rapid initial "saturation" mechanism is observed 
at an applied field of 400 Oe only. In the canting mech- 
anism discussed above, we have ascribed this low-field 
process with an alignment of all the spins on the individ- 
ual sublatticcs. This implies the removal by the field of 
all domain-walls along the chains, as well as the break- 
ing up of the interchain interactions. It follows that the 
effective field H' = z'J'S/gfiB should be smaller than 
400 Oe, which yields (taking g ~ 2, and S = 1) the up- 
per limit J'/fce ^ 0.03 K. This would suggest in fact that 
the breaking up of the interchain couplings occurs in even 
much smaller fields already and that the magnetization 
process up to 400 Oe mainly involves the removal of the 
domain-walls from the chains by the Zeeman energy. 

The quite peculiar form of the hysteresis in Fig. 10 
appears to agree with the just sketched scenario. Af- 
ter first saturating the chain system, the ferromagnetic 
order within and between the chains is seen to remain in- 
tact until the applied field is reduced to effectively zero, 
whereafter the magnetization discontinuously drops to 
zero. The instability of the ferromagnetic order in zero- 
field even at 1.8 K indicates that the interchain interac- 
tion is most probably antifcrromagnetic. As soon as the 
applied field has become smaller than this antiferromag- 
netic coupling, it will tend to order the ferromagnetic 
chains back again into an antiparallel arrangement. This 
means that half of the chains will have to switch their 
magnetization and as we have argued in the above, this 
can only occur by the excitation and subsequent propa- 
gation of domain-walls. Part of these walls will then be- 
come pinned again so that the antifcrromagnetic arrange- 
ment can not become complete and when subsequent by 
increasing the field in the opposite direction, the same 
cycle of removal of the walls has to be overcome. 
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VI. CONCLUDING REMARKS 

Clearly at variance with all previously reported behav- 
ior of both polytypes of iron(II) phthalocyanine, we have 
found that the magnetic structure of a-FePc can be de- 
scribed in terms of ferromagnetic Ising chains, with very 
weak antifcrromagnetic interchain coupling. As a con- 
sequence of impurities and other possible lattice defects, 
the transition to 3D ordering is suppressed, so that in- 
stead, a gradual freezing/blocking transition is observed 
to a disordered magnetic state. The disorder is lifted 
by a small field of about 2 kOe. Below 5 K this tran- 
sition is hindered by strong relaxation effects, which we 
attribute to the pinning of the domain-walls by the ran- 
domly distributed defects in combination with the in- 
terchain coupling. The competition between magnetic 
exchange interaction and crystal field effect gives rise to 
a canted arrangement of the ferromagnetically aligned 
spins. Above 10 K the behavior of the susceptibility 
can be interpreted in terms of the thermal excitation 
of kinks (domain-walls) and kink-antikink pairs. Pre- 
liminary Mdssbauer experiments can be consistently ex- 
plained in the same terms, the kink dynamics being re- 
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ent temperatures. Above 25 K, the shape of the 
spectra is almost independent of temperature. 

12. A schematic of a g -E g interaction for a- form (a) 
and /3-form (b), and e g -E g interaction for a- form 
(c) and /3-form (d). 
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